function y=Theorem2(u,mu,sigma,a,alpha)
y=(2*mu/sigma^2-1)*a^(2*mu/sigma^2-1)*alpha^(-2*mu/sigma^2);

